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Abstract: The increase in computation and sensing capabilities as well as in battery duration 
of commercially available Wireless Sensors Network (WSN) nodes are making the paradigm 
of an horizontal ambient intelligence infrastructure feasible. Accordingly, the sensing, 
computing and communicating infrastructure is set with a programmable middleware that 
allows for quickly deploying different applications running on top of it so as to follow the 
changing ambient needs. In this scenario, we face the problem of setting up the desired 
application in complex scenarios with hundreds of nodes, which consists of identifying 
which actions should be performed by each of the nodes so as to satisfy the ambient needs 
while minimizing the application impact on the infrastructure battery lifetime. Accordingly, 
we approach the problem by considering every possible decomposition of the application's 
sensing and computing operations into tasks to be assigned to each infrastructure component. 
The contribution of energy consumption due to the performance of each task is then 
considered to compute a cost function, allowing us to evaluate the viability of each 
deployment solution. Simulation results show that our framework results in considerable 
energy conservation with respect to sink-oriented or cluster-oriented deployment approaches, 
particularly for networks with high node densities, non-uniform energy consumption and 
initial energy, and complex actions. 
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Energy Consumption 



P T Transmitting radio frequency power consumption 

P R Receiving radio frequency power consumption 

P T0 Transmitting circuit power consumption 

Pro Receiving circuit power consumption 

P A Power Amplifier power consumption 

P Tx Output power at the antenna 

P Rx Input power at the antenna 

5 Distance between transmitting and receiving nodes 

i] Drain efficiency of the Power Amplifier 

A Parameter determined by the characteristics of the antennas such as g 

and efficiency 

a Path loss exponent 

e Parameter proportional to P Rx and A 

e tx Transmission energy consumption 

R Data rate 

6t Energy consumed to send one bit 

cr Energy consumed to receive one bit 

e sens Sensing energy consumption 

e proc Processing energy consumption 

e mstr Average energy consumption per instruction executed 

M task Number of instruction to process task 

smp m Number of samples to be processed 



Problem Formulation 



X = {xi, ...,Xi, xn} Set of N nodes in the WSN 

A = (ay) G R NxN Adjacency matrix 

A = (<%) G R NxN Distances matrix 

E = (eij) G ~SL NxN Matrix of the parameters e 

V% — {Pmi-, PT0i,Vt} Set of characteristic parameters of the node Xj 

F — {/ 0 , /u,, fw} Set of W tasks that can be performed by any node in the network 

Di = {da, d im , d^} Set of /j tasks that the node is able to perform 

O Operation performed by the WSN 

E tot Total cost value associated to O 

C — {ci, q, c L } Set of L sub-operations which must be executed by the nodes to perform 

operation O 

S = {si, sn} Set of N statuses, where Sj corresponds to the status of the node Xi 
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Traffic Flows 



4-OUt j h.out fout\ 

h ih — XHh i Jih i 



Vector of the output traffic flows generated by node Xi 

Traffic flow, corresponding to the h-ih. element of T° ut , where each 

sample of bits is transmitted at the frequency ffif 

Vector of the input traffic flows, corresponding to the union of the 

traffic flows received by node 

Function which generates the output traffic T° ut , based on the input 
traffic and the status of the node x, 



Virtual Nodes 



Gi 
L 

X vir = ^ x v 

A mr e I 

A vir e 

e I 



~\rvir 

v 



d: 



^ N vir xN vir 



rpVir \ 



Set of sub-operations that a single node Xi can perform 

consecutively for the implementation of the operation O 

Set of sub-operations that the nodes cannot execute with the others 

Set of N vir nodes of the WSN with virtual nodes 

Adjacency matrix for the WSN with virtual nodes 

Distances matrix for the WSN with virtual nodes 

Matrix of the parameters e for the WSN with virtual nodes 

Set of characteristic parameters of the node x v v ir of the WSN with 

virtual nodes 

Set of tasks that the node x™' of the WSN with virtual nodes is able 
to perform 



Cost Functions 

li 



j^proc 



QT 
QT 



Sensing cost function for node Xi 

Coefficient in inverse proportion with the residual energy of the 
node x,i 

Processing cost function for the node Xi 
Transmission cost function for the node x^ 

Minimum number of nodes that have to perform the sub-operation q 
Maximum number of nodes that have to perform the sub-operation q 



Introduction 



Over the past decade, Wireless Sensor Networks (WSN) have consistently evolved into more complex 
distributed monitoring and control systems. In the beginning, the only goal of a WSN was to monitor a 
given environment. As depicted in Figure 1, sensors gathered the required information mostly according 
to a fixed temporal schedule and sent it to the routers, whose primary functions were to both perform 
the sensing and conveniently transmit data received from sensors. Finally, the Coordinator would collect 
the data transmitted by the router. The Coordinator has much higher processing and memory capacity 
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than the other nodes in the network. It initializes and manages the WSN: it is responsible for routing and 
storing the routing tables as well as information about the network and security keys. The Coordinator 
usually interfaces with a server or a computer. In this way, data from sensors could be processed and 
stored. The Coordinator possessed all the intelligence of the network. Given the cost of more advanced 
devices, all other nodes only had basic processing and memory capacity [1]. 

A WSN is characterized by the presence of one or more sinks. A sink is a node which gathers and 
controls data collected by different sensor nodes. In every network there is at least one sink present, 
corresponding to the Coordinator. 

Figure 1. Structure of a WSN. 
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WSNs are now becoming more and more complex. According to the information gathered by 
the sensors, the network is capable of making decisions and acting upon them. Indeed, they are 
expected to be one of the pillars of the Internet of Things (IoT) paradigm [2], which foster the 
introduction of key applications, including but not limited to domotics, assisted living, e-health, 
enhanced learning automation and industrial manufacturing logistics, business/process management, and 
intelligent transportation of people and goods. Reduction in the cost of the devices has increased the 
nodes' capacity, thus they can perform some processing before sending the data to a sink. They can 
aggregate data coming from different sensors, perform temporal and spatial averaging as well as data 
filtering so as to reduce the burden of transmitting large amounts of data to the Coordinator and increase 
the network lifetime. Indeed, devices in a WSN are typically powered by batteries that can be difficult 
to replace, such as in the case of subterranean or underwater nodes. 

These considerations contribute to the vision of an horizontal ambient intelligence infrastructure 
wherein the sensing, computing and communicating infrastructure is set with a programmable 
middleware that allows for quickly deploying different applications running on top of it so as to follow 
the changing ambient needs, for example: monitoring a given geographical area and alerting when 
something is happening therein; activating the heating system when the ambient is getting cold; tracking 
the processing chain in industrial plants to prevent hazardous scenarios. In this case, we focus on the need 
for a logic that, starting from the desired application, can be set up in complex scenarios with hundreds 
of nodes, evaluate the possible deployment solutions and decide which action should be performed 
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by each of the nodes so as to satisfy the ambient needs while minimizing the application's impact 
on the infrastructure's battery lifetime. Accordingly, we approach the problem by considering every 
possible decomposition of the application's sensing and computing operations into tasks to be assigned 
to each infrastructure component. The contribution of energy consumption due to the performance of 
each task is then considered to compute a cost function allowing us to evaluate the viability of each 
deployment solution. 

This paper is organized as follows. The second Section describes how other authors have dealt with 
the problem of maximizing of the network lifetime, and the model of energy consumption considered 
in this work. The third Section introduces the problem and how we have approached it. The following 
Section defines the algorithm used by the framework to maximize the network lifetime. In the fifth 
Section, some simulation results on the effectiveness of the framework are presented. The last Section 
draws final conclusions. Finally, in the Appendix we provide the list of symbols used to make it easier 
to follow the description of the proposed algorithm. 

2. Background 

2.1. Past Studies 

Due to their scarce resources, minimization of energy consumption has been a key challenge for 
Wireless Sensor Networks. There are a great number of works that have focused on the maximization of 
the network lifetime, each taking into account a different approach to achieve it. 

Routing is probably the most immediate issue addressed in order to accomplish this goal. A 
convenient choice of paths to route data may result in significant energy conservation. In [3], an 
energy-efficient metric for finding routes was proposed. In [4], nodes energy reserves were taken into 
account to route the traffic so that the nodes' drain-out is maximized. Other routing techniques are shown 
in many other studies, such as [5] and [6]. 

Some studies build on the assumption that transmission energy consumption is related to the square 
of the distance between two communicating nodes. Therefore, it might be more energy efficient to send 
data over many short hops rather than fewer long hops. This issue was handled in [7] and, later, in [8]. 
This approach intends to maximize the network lifetime by minimizing overall energy consumption. 
However, it does not resolve the problem of unbalanced energy consumption among the nodes. This 
uneven energy dissipation may lead to an early death of some nodes, indeed resulting in a reduction of 
the network lifetime rather than an increase. 

In order for some nodes not to die much earlier than others, energy consumption in the network should 
be as balanced as possible. Relay nodes might be used for this purpose, as shown in [9] and [10]. Energy 
load distribution can also be achieved by conveniently deploying the network nodes, as in [1 1]. In [12], 
the nodes are spaced non-uniformly as a function of their distance. Taking into account that nodes near 
the base station feel the effects of higher traffic more than other nodes, spacing is adjusted in such a way 
that nodes with higher traffic have a shorter hop distance than nodes with less traffic. 

None of the studies mentioned above considers the possibility of processing the data in the nodes of 
the path to the destination. Because most of the energy spent in a Wireless Sensor Network depends 
on the amount of data that is transmitted over the network, reducing the amount of data may result in 
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a reduction of the transmission energy consumption. This principle has been only partially adopted 
by LEACH [13], where sensors serve as Cluster Heads aggregating the data and, indeed, decreasing 
the number of bytes sent over the network. Energy consumption balancing is guaranteed by a random 
rotation of the role of Cluster Head. 

In [14], Cluster Heads for a mobile WSN are chosen according to a selection algorithm which 
takes into account characteristics of the nodes such as power energy and transmission rate. The 
clustering algorithm presented is periodically re-run to possibly find new nodes which might improve 
the network lifetime. 

Given the computational capacity of modern sensors, a step forward could be taken not just by 
aggregating data, but by processing them before they arrive at their destination whenever possible and 
on the basis of the network topology and power resource detection. In this paper, a framework that 
determines which of the nodes should process the data in order to maximize the network lifetime is 
presented. To the best of the authors' knowledge, no similar frameworks for WSN have been proposed 
before. An example of an overlaying framework that handles an architecture for an integration of the 
IoT in enterprise services might be found in [15]. However, this framework is not conceived to minimize 
the energy consumption, unlike the framework hereafter described. 

2.2. Energy Consumption 

Energy consumption in Wireless Sensor Networks is determined most of all by transmission and 
reception. As mentioned in [16] 

Pt(S) = Pro + Pa(S) = Pro + PrMh (1) 
Pr = Pro 

where P T and Pr are radio frequency power consumptions for transmitting and receiving respectively; 
Pa is the power consumption of the Power Amplifier (PA); 5 is the distance between the transmitter and 
the receiver; P T0 and P R0 are the components of power consumption of the transmitting and receiving 
circuitry respectively; Pt x is the output power at the antenna which, for reliable transmissions, depends 
on the distance 5; r\ is the drain efficiency of the PA. 

Considering a channel in which the path loss component is predominant, and thus secondary effects 
such as multipath and Doppler can be neglected, the received power P Rx can be expressed as 

Prx = (2) 

where A is a parameter determined by the characteristics of the antennas (such as gain and efficiency) 
and a denotes the path-loss exponent, which is about 2 for free space. This kind of modelling is typical 
of free space propagation. Of course, the model might be extended to account for other fading effects. 
From Equations (1) and (2) 

/ c \ „ Prx x A x 5 a 
P T {5) = P T0 + 

V 

Considering e = P Rr x A, where P Rr is the minimum reception power for a 
reliable communication 

£ X S a 

Pt(S) = Pro + 

V 
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This implies that the total power consumption for communicating between a transmitting node A and 
a receiving node B of a WSN could be written as 

P - P ip i £ab x S ab 
±AB — -TROs + -TTOa "i 

Therefore, the energy consumption of the network to transmit a packet of k bits from A to B with a 
constant data rate R is 

txn r> d x \ -Pas x k 
e [K, ^T0 A ,VA, ^R0 b i £ AB,0ab) ~ ^ 

_ k ( e AB x 5% \ 

= k x (e RB +e TA (S)) (3) 

where e TA is the energy to send one bit over a distance 5ab and e# s is the energy to receive one bit. 

The model described does not take into account mechanisms such as sleep schedule and route 
discovery, which may produce overhead. Thus it could thus be necessary to consider not just the single 
packet transmission, but also the energy consumption due to the overhead. 

Besides transmission and reception, the other two main causes of energy consumption are due to 
the sensing activity and to the processing. The sensing energy consumption e sens is determined by the 
specific characteristics of the sensor, which is given by the used hardware. 

The processing energy consumption e proc (task, data m ) is proportional to the complexity of the 
task — that is, the number of instructions needed to complete the task — and to the ingress data 
data in — the higher the number of samples involved in the processing, the higher the energy consumption. 
This function e proc (.) can be determined on the basis of the device datasheet used. Calling M task the 
number of instructions for the task, smp m the number of samples to be processed and e mstr the average 
energy consumption per instruction executed 

e proc {task, data 111 ) = smp m x M task x e instr (4) 

A summary of the notation used throughout this paper is provided at the beginning of this document 
for quick reference. 

3. Problem Formulation 

The goal of a WSN is to accomplish a given number of operations mostly based on some 
measurements performed on the relevant environment. In our scenario not all the nodes have the same 
capacities, as represented in Figure 2, where three sets of possible tasks have been considered, which 
could be: data processing, temperature measurement and video monitoring. Given the status of the 
network in terms of node capacities, topology, and energy distribution, the problem addressed is to 
assign to each node the tasks that, combined together, contribute to the target network operations while 
minimizing the application impact on the infrastructure battery lifetime. 

In our modelling, the number of nodes X = {x±, Xi, x^} in the WSN is denoted with N, where 
the node Xi can be a sensing node, a router or an actuator (or node with a combination of these roles). 
The node x N refers to the sink (we assume there is only one sink in the network). The network can be 
described by: 
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• the N x N adjacency matrix A = (ay): an element of A is equal to 1 when a link connects 
node i to node j and the sink is closer to j than to i; 

• the distance matrix A = (5^), which contains the pairwise distances (in meters) between adjacent 
nodes. If Sij = oo, nodes i and j are not adjacent; 

• the matrix E = (%•), with the parameters introduced in Section 2.2, calculated for each 
combination of adjacent nodes i and j, where i reaches the sink though j. If = oo, nodes 
i and j are not adjacent; 

• the set of characteristic parameters Vj, = {Proi, Proi,Vi} °f every node x i5 which are useful to 
compute the transmission energy consumption as defined in Equation (3) in Section 2.2; 

• the set F = {f 0 , fw} of tasks, which encompasses all the tasks that can be performed 
by any node in the network. For instance, 0 might correspond to "temperature sensing in the 
area 1", 1 to the "temperature sensing in the area 2", 2 to the "pressure sensing in the area 3", 
3 to the "spatial averaging" (which means performing the average of the sensed data arriving 
from different geographical areas), 4 to the 'temporal averaging" (which means performing the 
average of sample values sensed by the same sensor at different instants of time), 5 to the "only 
transmission", 6 to "no actions". Each of the tasks in an F set entails a transmission of data, with 
the exception, of course, of "no actions"; 

• the set Di = {da, d im , du.}, with Di C F, where the elements of Di are the tasks that the 
node Xi is able to perform. 

Figure 2. Example of a WSN. Nodes belonging to set 1 might perform task fi, nodes in 
set 2 might perform task f 2 \ nodes in set 3 might perform task f 3 . 




We assume that a given operation O has to be deployed in the network, which can be decomposed 
into a sequence of distributed tasks. This could represent diverse operation, including: computing 
the average of the temperature in certain geographical areas, measuring the light intensity in a room, 
video-surveillance of a specific geographical area, or a combination of these. In the following, we rely 
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on a specific reference application, which we name spatial and temporal monitoring: a spatial and 
temporal mean operation over an hour is performed on the temperature values sensed by the sensors 
every 10 min by the sensors from 3 different locations; the average values are stored in the sink. We use 
this example to better explain our modelling. 

Three significant parameters can be associated with the operation to be deployed: 

• the total cost value E tot , which takes into account the energy consumption; 

• the sequence C = {ci, q, Cl] of the sub-operations that must be executed by the nodes to 
perform the operation O. If, for instance, O is our reference example of spatial and temporal 
monitoring, some of the sub-operations in C will certainly be "temporal averaging" and "spatial 
averaging". Each sub-operation is equivalent to one of the f w tasks defined above: C C F. 
The terms sub-operation and task are used interchangeably, but we refer to q as a sub-operation 
to emphasise its link with its respective operation O. Sub-operations are listed in C in order of 
priority: if a node has to execute both c\ and c 2 , the former must be executed before the latter. 
With reference to our spatial and temporal monitoring example, c\, c 2 and c 3 are the sensing 
sub-operations: c\ is "temperature sensing in area 1", c 2 is "temperature sensing in area 2" and C3 
is "temperature sensing in area 3"; because every sensor only belongs to one area, the order of the 
sensing sub-operations is irrelevant. Rather, it is important that data is gathered before any other 
sub-operation is performed on it; thus, the sensing sub-operations are the first elements in C. If 
a node has to calculate both the temporal and spatial mean of the received values, it has to first 
work out the temporal mean on the data received from each path, and then it has to perform the 
spatial mean on it. Therefore, C4 is "temporal mean" and c 5 is "spatial mean". Finally, c$ is "only 
transmission" and C7 is "no actions". All the sub-operations for the example are summarized in 
Table 1; 

• the set S = {s 1 , Sj, sn}, where Sj is the status of node i with respect to operation O. The 
value Si defines which sub-operation q the node i is performing. Of course, the status of the node 
%i has to be chosen among the set of tasks Di that the node is able to perform. If the node is not 
involved in the operation O, its status corresponds to "no actions". For this reason, "no actions" is 
necessarily included in set F. The following must always be verified: 

- if a node Xj receives some data from a node Xi, which means that Sj cannot be "no actions", 
node Xj must at least transmit the data 



- if a node xj does not receive any data, which means that the status of all the nodes connected 
to it is "no actions", Sj must also be set to "no actions" 



Si 7^ "no actions" V 



1 =>- Sj ^ "no actions' 



(5) 



Si = "no actions" Vi : a,ij = 1 Sj = "no actions 



(6) 
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Table 1. Sub-operations for the spatial and temporal monitoring example. 





Flpcf i*SntS (\w 


Cl 


Temperature sensing in area 1 


C2 


Temperature sensing in area 2 


C3 


Temperature sensing in area 3 


c 4 


Temporal mean 


c 5 


Spatial mean 


C6 


Only transmission 


C7 


No actions 



Thanks to the greater processing power and storage capacity of modern sensors, contrary to the past, 
the same operation O can be performed in several different ways: gathered data can be immediately sent 
to a sink or it can be processed before being transmitted. In the case of the latter, the number of bits to 
be sent would be smaller, and therefore the transmission energy consumption would be lower as well; 
however, processing energy consumption could be higher in this second case. Quantifying the energy 
consumption in both cases, it could be possible to establish which one determines a reduction of battery 
consumption in the sensors, incrementing the network lifetime. 

The framework described further takes a high-level code as input, evaluates which combination of 
statuses S = {si, s i: s N } permits the operation O to be performed with the lowest possible energy 
consumption E tot , and finally elaborates and assigns among the nodes the most appropriate tasks to be 
performed. Hence, it is evident that the cost function E tot will vary depending on the status of each node, 
that is, how the operation O is performed. The problem addressed is then defined as the set of statuses S 
that minimizes the impact of the application on the network lifetime. In the following, we elaborate the 
considered scenario by defining further constraints that solve the problem. 

4. Deployment of Distributed Applications 

In the following, we present the proposed solution towards a distributed application deployment in 
WSN. The following Subsections present: the constraints on the traffic generated by the distributed 
applications; the concept of virtual nodes, which are duplicates of real nodes that are introduced to deal 
with nodes that perform more than a single action in a deployment solution; the cost functions built on 
the basis of the energy consumption formulae; the network lifetime maximization procedure; a summary 
of the proposed framework. Note that the modelling that we propose in this work and that we present in 
this Section is aimed at evaluating all the possible solutions of application deployment in terms of data 
transmission and processing. The parameters, constraints and cost functions are introduced for the sole 
aim of evaluating the viability of the solution. 

4.1. Constraints on Traffic Flows 

In our scenario we assume that the sources of traffic in the network (the sensors) generate samples of 
k bits at a certain frequency /. The processing in the network is performed on this type of traffic flow 
coming from different nodes. The generic node xi receives the traffic T* n over which it performs the 
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sub-operation corresponding to its assigned status s ?; . The effect of this sub-operation is the generation 
of the output traffic T° ut , which is computed by function p as follows 

Tr*=p(Tf, Si ) (7) 

The output traffic is then sent to the next node towards the sink. 

The data generated by p in the node is modelled by the iJ-dimensional vector 
T° ut = (t«* £%*,..., t$? ), where each element t$* = {k$*, /X*} corresponds to a traffic flow where 
each sample of k°^ 1 bits is transmitted at the frequency f™*. Each sample described by tffi results from 
a spatial processing or a sensing. The data T° ut is then sent to the following node Xj, according to 
adjacency matrix A. 

The node Xj receives data from all adjacent nodes that reach the sink through Xj, with the exception 
of the nodes with a "no actions" status 

N 

TJ 1 = {jT° ut xdijx z u with Zi = 
i=i 

As defined by Equation (7), the data TJ 1 received by the node Xj is processed, according to the status 

of Xf 

• if Sj is a sensing status, p does not take any TJ 1 as input and the output is defined by the specific 
sensing operation; 

• if Sj is an "only transmission" status, the output of p is exactly equal to TJ 1 ; 

• if Sj is a "no action" status, p returns a T° ut with all fields set to 0 as output. This case is included 
for completeness, but it is not supposed to happen because of Equation (5); 

• if Sj is a processing status, T° u * can be the most diverse depending on the specific processing 
objectives, which are coded in Sj and that control the specific function p. In the following we 
analyze certain cases. 

Referring to the spatial and temporal monitoring example, processing can be a spatial averaging, a 
temporal averaging, or a combination of both. In a spatial processing, the samples coming from different 
paths are processed together, as shown in an example in Figure 3(a). Here, four flows of 25 bits per 
second are received by node 4, which are then averaged to produce a single flow of 25 bits per second. 
Accordingly, the resulting T° ut is made of only one element t™ 1 = {k^ f , where the number of 

bits per sample k°^ and the frequency are equal to those of each input flows. It must be noted that, 
in general, k°^ is not necessarily equal to the number of bits of each input flow, but it may be different 
according to the processing output. 

Differently, the temporal averaging is performed on every traffic flow in T* n . The resulting T° Mi 
contains the same number of traffic flows as in TJ 1 , where every element tfffi is characterized by the same 
number of bits per sample k°~jf and the same frequency corresponding to the averaging frequency 
associated to the node status Sj. Indeed, we may have different status codes associated to the temporal 
averaging, each one distinguished by a different processing frequency. Figure 3(b) shows an example for 



0 Sj = "no actions" 

(8) 

1 otherwise 
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this kind of processing, where, in this case, we assume the node status corresponds to temporal averaging 
with frequency 0.5 Hz. 

Other processing tasks can be performed on every single sample of each received traffic flow without 
involving other samples. This is the case, for instance, where one must evaluate whether the received 
values exceed a given threshold or not, consequently transmitting a boolean output value. The only thing 
that changes in the output traffic flows is the number of bits per sample; therefore, T° ut contains the 
same number of traffic flows as in T* n at the same frequency fjjf, but with different bits per sample kpf. 
Figure 3(c) shows the traffic flows for the described processing. 

Figure 3. Examples of data processing in node 4 that receives input traffic from nodes 1-3. 
In these sketches we show the input and output traffic for: spatial averaging (a); temporal 
averaging (b); and single sample processing (c). 




T° 2 ut = [25,2] 



(c) Single sample processing 

There are many other processing tasks that can be performed in a given network and they are coded 
in F as described in the previous Section. For each one of these, an operator p(x, y) is defined. Note that 
for our objective, this operator is needed to figure out the traffic flows that will be traversing the network 
for each deployment scenario. 

4.2. Virtual Nodes 

It is possible that a single node has to perform more than one task. For instance, if operation O is a 
temporal and spatial average of the temperature values measured in different geographical areas, it may 
happen that a single node has to compute both spatial and temporal average values on the received data. 
To take into account this type of scenario, we rely on the concept of virtual nodes. These are copies 
of real nodes, each one able to perform only a specific sub-operation and sending the resulting data at 
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zero-energy cost to the next virtual node (except the last one that sends the data to the next node). We 
define the set of sub-operations that a single node Xi can perform consecutively for the implementation 
of the operation O 

Gi = Di f] C \ L (9) 

where L is the set of sub-operations that the node cannot execute with the others, which are "only 
transmission" and "no actions". These two are kept outside the set Gi because it cannot happen that a 
sequence of sensing/processing tasks are followed by any of the tasks in L. If > 1, the node Xi is 
divided into \Gi\ virtual nodes. Each virtual node is created so as to be able to execute only one of the 
sub-operations in Gi, to which the sub-operations in L are added. Each virtual node can then be assigned 
to perform only one of these three sub-operations (the one taken from Gi plus the two in L). Hence, the 
new network will have a number N mr of total nodes x™ 

N 

N vir = Y^ Tl , Withr i = 

i=l 

The set of possible sub-operations for x v v tr is D^ ir . Figure 4 draws an example of sequence of virtual 
nodes for node x±, which is substituted by nodes x\ ir , x™ r and x v § T . Let us refer to the spatial and 
temporal monitoring example and the associated set of sub-operations C listed in Table 1 . Additionally, 
let us assume that {ci, c 4 , c 5 , c 6 , c 7 } C D 4 , that is: the node x 4 can sense the area 1, perform a temporal 
averaging, a spatial averaging, only transmit the received data and perform no actions (these last two 
tasks are always included). Because L = {cq, c 7 } 

DZ = {c u c 6 ,c 7 } 

Df r = {c 4 ,c 6 ,c 7 } (10) 
Dr = {c 5 ,c 6 ,c 7 } 

A new adjacency matrix A mr is defined to incorporate additional virtual nodes. Such matrix is built 
simply by substituting the real node with the sequence of virtual nodes, so that the first virtual node is 
connected to the nodes from which Xi received the data, while the last virtual node is connected to the 
node to which X{ sent the data. The other nodes are connected in sequence. An exception happens if the 
real node can also perform some sensing functions. In this case, the corresponding virtual node is kept 
outside this sequence and it merely sends the data to the subsequent virtual node. This rule has been 
introduced because the sensing operation does not need any data from other nodes. This scenario is what 
happens in the example shown in Figure 4, where the temperature sensing function can be performed by 
node x™ r . 

As defined in Section 3, sub-operations are in priority order. For this reason, virtual nodes have to 
be placed such that the first performs the sub-operation with the highest priority, and the last executes 
the sub-operation with the lowest priority and sends the data to the rest of the network. With reference 
to Figure 4, the highest priority is given to the sensing action, because whether there are any data to 
be sensed, the information should be gathered before processing it. For this reason, sub-operation c\ 
is included in D v " . If data received by node x™ r has to be processed, temporal averaging should be 
performed before spatial averaging. Thus, priority of the former must be higher than the latter, and this 



if|Gi|>0 
otherwise 
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is the reason why temporal averaging c 4 is an element of set D^ ir , while spatial averaging c 5 is an element 
of set Dq it . The adjacency matrix A mr has to be defined accordingly. 

Figure 4. Example of virtual nodes substituting node x±. Node x^ can perform a 
sensing action. 




x v 3 ir Ex 3 



A set of characteristic parameters V™ r = {P^, Ptov nas to be associated with each node. 
Taking into account the virtual nodes substituting node Xi, the transmission cost from a virtual node to 
the other must be null. Therefore, only the last virtual node, that has to transmit the data to the network, 
has the same characteristic parameters as x^. the virtual nodes before it must have the parameters 
and P™ set to 0. 

New matrices A w and E w are defined, where the values of their elements for adjacent virtual nodes 
from the same original node are null. From now on, we will be referring only to a network with virtual 
nodes; however, to make the presentation clearer, we will skip subscript "t>zr", as it is unnecessary. 

4.3. Cost Functions 

The objective of the proposed algorithm is to evaluate the viability of each deployment solution on 
the basis of a cost function that is connected to energy consumption. Quite often in similar scenarios, 
past studies have proposed the evaluation of the network lifetime and have aimed at maximizing it. 
Since in our framework we assume the network may be employed to perform more than one operation 
simultaneously, there is no sense in computing the network lifetime since it is affected by other 
applications which are not considered in the same analysis. For this reason we try to minimize the 
energy consumption for the operation under analysis, allowing the network administrator to also include 
a parameter that takes into account the current node battery energy level, as it is shown in the following. 

We consider three cost functions: one for the sensing, one for the processing and one for the 
transmission. The sensing cost function for node Xi is expressed as 

1 if Si = sensing code 

E sens = f out x ^ . x ^ens x ^ ^ y . = J 6 (n) 

I 0 otherwise 

with e| ens representing the sensing energy consumption as defined in Section 2.2. Recall that f° ut is 
the node output traffic frequency, which also represents the sensing frequency. The parameter 7, is 
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a coefficient in inverse proportion to the residual energy of the node which can be set to drive the 
deployment of the application towards nodes with higher residual energy levels, as anticipated above. 
When performing our experiments, we have set 7, to 1 when the battery is fully loaded, while we set 7$ 
to 5 when the battery level is lower than 20% of the total charge. From 1 to 5, 7, changes linearly. 
We defined the processing cost function as follows 



h= 



otherwise 



where 7$ is the coefficient defined above, and ef™ c is the processing energy consumption defined by 
Equation (4), which depends on the characteristics of the node Xj, the sub-operation c Si that has to be 
executed, which in turn depends on the status Sj of the node, and the received data T* n described in 
Equation (8). Because the processing cost depends on the number of processing per second performed 
by the same node Xi, it is proportional to the frequency f°£ l of each of the H egress traffic flows, where 
H is the size of T° ut as described in Section 4. 1 . The number of samples to calculate ef roc is defined 
differently for each kind of processing detected in Section 4.1. For a spatial processing, the number 
of processed samples is equal to the number of ingress traffic flows; for a temporal processing, the 
number of processed samples for each traffic flow is the ratio between the frequency of arrival of the 
samples f-fi and the processing frequency /^**; by definition, for a single sample processing the number 
of processed samples is 1. 

Both sensing and processing are followed by a transmission. Therefore, unless the node is in a "no 
actions" status, it has to transmit the output data. Because for Equation (5) if a node receives data it 
cannot be in "no actions" status, every node involved in the operation O has to transmit data. The related 
cost function is 

N N N 

Ef = fi x 7l x e tx (TT\ V h (J V, x a ij} (J e {j , (J c%) (13) 

j=i j=i j=i 

where is the transmission frequency, 7, is the residual energy coefficient, and e tx is the transmission 
energy consumption defined by Equation (3), depending on: the data to be transmitted T° ut ; the 
characteristic parameters Vi of the node xf, the characteristic parameters Vj of all the Xj nodes that 
will receive the data from x^ which, for a connected graph, is just one; the parameter £„• concerning 
nodes x ri and xf, the distance 5^ between xi and xj. 

Given Equations (1 1-13), the overall cost function for any operation O is 



iV 



E tot = ^ (E* ens + Ef roc + E\ x ) (14) 



8=1 



4.4. Maximization of Network Lifetime 



The goal is to find the set of the statuses S = {si, Sf, sjy} of the nodes that minimize the 
network energy cost function. 
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Therefore, the optimization problem becomes 



N 



minimize E tot = ^ (E? e ns + Ef oc + Ef ) (15a) 



i=i 

N 



subject to QT in <^qa < QT ax with ^ = J 1 Si ° l (15b) 
i=1 I 0 otherwise 

\^J{s v } = Gn \/v:x v ^x N (15c) 

V 

The condition in Equation (15b) is a constraint on the minimum (Qf im ) and the maximum {Q™ ax ) 
number of nodes that have to perform the sub-operation q. This could be necessary, for example, when 
a given geographical area is monitored by a certain number of nodes, but the required information is not 
needed from all of them. If, for instance, the mean temperature value of an area monitored by 30 sensors 
is needed, we may prefer the temperature information to be gathered just by 10 of those sensors, in order 
to consume less energy. In this case, both Q™ ln and Qf wx would be equal to 10, and the algorithm would 
choose the 10 sensors which weight less on the network lifetime, among the 30 sensors which are able 
to sense temperature in the required area. When this constraint is not needed for a sub-operation /, Q™ m 
is null and Q™ ax is set to N. 

The condition in Equation (15c) shows that the set of statuses of the virtual nodes corresponding to 
the original sink node x^ must correspond to the set |Gjv|- 

This implies that none of the virtual nodes corresponding to the original sink node can be in "only 
transmission" or "no actions" status, but they are in a processing status; therefore, if there is any data 
still to be processed, those virtual nodes have to process them. 

The problem defined in Equation (15) is a Binary Integer Programming (BIP) problem: the unknown 
status of a node X{ can be defined as a \C\ -dimensional binary array, where C is the set of sub-operations 
as defined in Section 3. Because every node can only have one status, which means that it can perform 
only one sub-operation among those that it is able to perform, only one element of this array can be equal 
to 1, and it corresponds to one of the sub-operations that the relating node Xi is able to perform, according 
to A- The elements of the array represent the weights to the contributions in Equations (11,12) of the 
node to the cost function. 

BIP problems are classified as NP-hard. Their exact solution is usually found using branch-and-bound 
algorithms. The worst case complexity of branch-and-bound algorithms is the same as the complexity 
of exhaustive search, which means that its complexity scales exponentially with the problem size. 
In the case under consideration, the problem size is dominated by the number of sub-operations \C\ 
and the number of virtual nodes N mr . Therefore, the worst case complexity would be 0(2^ xNmr ). 
Nevertheless, in most cases branch-and-bound is more efficient compared to exhaustive search. 
Furthermore, our problem's structure is such that only one element of the |C*| -dimensional array 
representing the status of each node is nonzero. This condition allows to reduce the search space to 
0(N mr ' c '). It has to be noted that commercial mathematical programming solvers such as CPLEX [17] 
or Xpress Optimization Suite [18] are claimed to use optimized branch-and-bound algorithms whose 
complexity scales linearly with the problem size. 
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In order to further reduce complexity, heuristic algorithms might be used as well, obtaining 
sub-optimal solutions which may be considered sufficient in most cases. 

4.5. The Proposed Framework 

Given a network similar to the one in Figure 2, the Coordinator, which is responsible for initiating and 
maintaining the network, is the device on which the deployment algorithm is performed. The proposed 
algorithm needs to know the exact topology of the network, that is, how the nodes are connected to each 
other and what the distance between any two of them is, as well as the routing table. In order to compute 
the cost function in Equation (14), further information is needed, such as the parameters to model the 
radio channel, the transmission, reception, sensing and processing energy consumption of each node, the 
residual energy of each node, the working frequency and the data rate. 

In summary, the algorithm performs the following steps: 

1. define iV Xj, sets F, Di, Vi and matrices A, E and A; 

2. define set O and the corresponding set C ; 

3. define sets Gi, 

4. define the new network with N vir virtual nodes xf r , and new sets Df r , Vf" ', A vir , E mr and A vir ; 

5. given in Equation (15), solve it with a linear programming solver, in order to find set S. 

The solving algorithm has been implemented in Mosel language, and the solution has been found 
using Xpress Optimization Suite. The binary array £™ r = {a™{, <r^ r , <J v v l [} has been associated 
to each node x v v ir ', where L is the cardinality of C, that is, the number of sub-operations which must 
be executed to perform operation O, as described in Section 3. The elements of E™ r must satisfy the 
following constraints: 

• the element a^f can be equal to 1 if and only if x v v tr is able to perform the sub-operation q e C, 
which means that cr^"' cannot be equal to 1 if the sub-operation q is not an element of the set of 
tasks D v v tr that the node x™ r is able to perform 

q i DZ aj ^ 1 (16) 

• only one element in can be equal to 1 

L 

Y,°vi = 1 ( 17 ) 

i=i 

The elements of the array built this way are the weights %)i and Vi of the energy contributions in 
Equations (11,12) defined in Section 4.3. 

The optimum way to perform the operation O at issue and spend the least amount of energy as possible 
will then be found. The node or the combination of nodes that are able to perform it and consume the 
minimum amount of energy will be chosen; then, a low level code describing which tasks each node has 
to perform will be developed and distributed to the appropriate nodes. 
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5. Performance Analysis 

5. 1. Test Cases and Simulations Setup 

To evaluate the effectiveness of the proposed strategy, three test cases have been taken into account 
according to some of the most significant realistic scenarios considered in past works, such as in [13]: 

1. easel: uniform energy consumption and uniform initial energy (UC-UE) at each node (equal 
characteristic parameters and battery life for every node); 

2. case2: non-uniform energy consumption and uniform initial energy (NUC-UE) at each node 
(different characteristic parameters but same battery life for every node). The energy consumption 
of the nodes have been assigned according to a uniform distribution from 60% to 140% of the 
energy consumption for case 1; 

3. case3: uniform energy consumption and non-uniform initial energy (UC-NUE) at each node (same 
characteristic parameters but different battery life for every node). The battery charge has been 
assigned randomly according to a uniform distribution from 20% to 100% of the total charge. 

We monitored a rectangular-shaped outdoor environment (e.g., a vineyard, a seaport, a tourist plaza), 
divided into areas of 25 m 2 , where the nodes have been deployed with different densities: 

• 0.2 nodes per square meter; 

• 0.3 nodes per square meter; 

• 0.4 nodes per square meter. 

The nodes have been positioned randomly, following a uniform distribution. Each node is equipped 
with sensors for the measurement of temperature, humidity, PH and light exposure. The data are sent to 
the Coordinator, which has identification number N. 

We have focused our analysis on the following two operations: 

(OpA) calculation and storage in the sink of the (temporal and spatial) mean values of temperature, 
humidity, PH and light exposure over an hour, starting from the values sensed every 10 min from 
every area; 

(OpB) aggregation of traffic coming from different areas of the network, carrying temperature, humidity, 
PH and light exposure values to the Coordinator for later analysis by qualified staff. 

We assumed that each sensed value is represented as a double numerical value, which is 64 bits long. 
Note that these two operations have been selected to compare the scenarios in which the network is 
required to perform significant data processing (OpA) and nodes have to perform only basic processing 
on the data, yet can significantly reduce the amount of transmitted data by aggregating the sensed 
samples (OpB). 

For both OpA and OpB, the first sub-operations are the sensing ones: c\ is the "temperature, humidity, 
PH and light exposure sensing for area 1", c 2 is the "temperature, humidity, PH and light exposure 
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sensing for area 2", and so on for every area. In order not to weigh down the text with alternatives, we 
suppose there are just 2 areas, which is not the case of the simulation scenario. Similarly to the spatial 
and temporal monitoring example in Section 3, for OpA c 3 is "temporal mean", c 4 is "spatial mean", c 5 
is "only transmission" and c 6 is "no actions". In addition to C\ and c 2 , the sub-operations for OpB are: 
C3 as "aggregation of samples", c 4 as "only transmission" and C5 as "no actions". The sub-operations for 
the two operations are summarized in Table 2. 

The nodes communicate using IEEE 802.15.4 radio interfaces on the 2.4 GHz frequency band. To 
keep things simple, any possible overhead has not been taken into account. 

We have simulated the resulting scenarios in MatLab environment, where we have implemented the 
proposed algorithm together with alternative approaches as discussed in the following Subsection. The 
main setup parameters are listed in Table 3 ( [16,19,20]). 



Table 2. Sub-operations for OpA and OpB, for two monitored areas. 





OpA 


OpB 


Cl 


Temperature, humidity, PH and light exposure sensing for area 1 


C2 


Temperature, humidity, PH and light exposure sensing for area 2 


C3 


Temporal mean 


Aggregation of samples 


C 4 


Spatial mean 


Only transmissions 


C5 


Only transmission 


No actions 


C6 


No actions 





Table 3. Simulation setup parameters. 



Parameter Value 



RF frequency 2,400 MHz 

Bit rate 250 kbps 

Programmable output power range Programmable in 8 steps from approximately —24 to 0 dBm 

Receiver sensitivity —94 dBm 

P m 59.2 mW 

P Ta 26.5 mW 

77 50% 

A 14 dB 

e instr 1 nJ 

Packets header 12 bytes 

Packets maximum payload 125 bytes 



5.2. Analysis of Case Studies 

The optimization algorithm has been applied to each of the cases mentioned in Section 5.1. The 
resulting cost value has been compared with: 

1. the cost value that is obtained if data are processed only by the Coordinator. This means that 
each traffic flow generated by the sensors is sent to the Coordinator without any processing at the 
intermediate nodes. We refer to this comparison with the letter "C"; 
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2. the data are processed (whenever needed) by every Cluster Head (node receiving flows from 
different sensors) found in the path to the Coordinator. We refer to this comparison with the 
letters "CH"; 

3. the mean cost value for all possible solutions that might be detected. This is introduced to make 
a comparison with a possible solution where the processing of the data is performed on fixed 
nodes, which is expected to bring results corresponding to the median solution. We refer to this 
comparison with the letter "M". 

These comparisons are expressed as a percentage of the energy conservation that would result using 
the proposed technique with respect to the alternatives one. 
Figure 5 shows the results for the two operations. 

Figure 5. Percentage of energy conservation using the framework, for (a) OpA and (b) OpB, 
in cases UC-UE, NUC-UE and UC-NUE, comparisons C, CH and M, for a node density 
of 0.3 nodes per m 2 . 




(a) OpA 




(b) OpB 



Sensors 2011, 11 



7415 



The two graphs show significant improvements of the proposed strategy with respect to the alternative 
ones with an average improvement of 36.3%. Limited benefits are observed in case UC-UE for both 
OpA and OpB. In fact, when all the nodes have the same parameter, and thus have a uniform energy 
consumption and the same initial energy, the choice of which node will perform the processing boils 
down to which Cluster Heads will do it. To illustrate the scenario, we refer to Figure 6, which depicts 
a Cluster Head CH 1 connected to the Cluster Head CH 2 by nodes 1, 2, and 3. Note that the Cluster 
Heads are just nodes that receive more than one traffic flow from different links. Because processing in 
CH 1 weights on the network as much as processing on node 1, node 2 or node 3, any processing of the 
data before arriving to CH 2 is more energy conserving. Processing the data on CH 1 ensures spending 
less transmission energy than processing data on nodes 1, 2 or 3. Accordingly, the CH approach allows 
for obtaining results similar to the ones we obtained with our approach in case UC-UE. We get slightly 
better results because when Cluster Heads are close to each other sometimes it is better to perform the 
processing only in the second Cluster Head rather than in both of them. 

Figure 6. Example of a transmission from Cluster Head CH 1 to Cluster Head CH 2. 




On the contrary, in cases NUC-UE and UC-NUE, devices' energy consumption does not weight the 
same amount on the entire network. This means that the nodes chosen by the proposed algorithm to 
perform the processing will be those that weight less on the network, regardless of whether they are 
Cluster Head or not. Therefore, the detection of the lower cost solution determines the best results, in 
terms of energy consumption, for networks with heterogeneous parameters, which are the most common 
type of networks in real scenarios. 

The benefits with respect to the approach CH are lower than the case of using approaches C and M. 
In fact, literature dictates that the use of Cluster Heads is a convenient solution, because the aggregation 
of frames coming from different paths leads to a reduction in network energy consumption. For this 
reason, when using Cluster Heads the cost is much lower, compared to sending every single frame to the 
Coordinator, or to the average of the other possible solutions; this determines less difference from the 
optimization algorithm solution, and thus a lower energy conservation. However, this approach requires 
every node in the network to be able to perform data processing, which is not always the case. In any 
case, our approach is proven to always outperform the CH approach. 

It could be noted that for OpA, in which processing is more elaborate and the number of instructions 
for every process greater, energy conservation is higher than it is for OpB. In fact, as could be expected, 
the lower the energy cost necessary for the processing, the more convenient it is to process the data in 
every Cluster Head encountered. 

This fact is demonstrated by the results shown in Figure 7, which depicts the percentage of energy 
conservation while the ratio between the processing cost and the cost to transmit 137 bytes of data 



Sensors 2011, 11 



7416 



increases. The distance considered is equal to the average distance of all the nodes from the Coordinator. 
Comparison has been made both in the case that data are processed only by the Coordinator (solid lines) 
and cases in which data are processed by every Cluster Head (dashed lines). In the former, energy 
conservation decreases when processing cost to transmission cost ratio increases. In fact, as could be 
expected, when the processing cost increases with respect to the transmission cost, it becomes more 
convenient to transmit data rather than process data. On the other hand, when compared with the CH 
approach, energy conservation increases when the processing cost to transmission cost ratio increases. 
In fact, when the processing cost increases, it is more convenient to accurately choose the nodes where 
processing might be performed rather than processing data every time that it is possible to do so. Table 4 
shows the results for OpA and OpB, for different node densities of 0.2 and 0.4 nodes/m 2 . The resulting 
tendency of an improved energy conservation when node density increases is basically due to two factors: 

• in cases NUC-CE and UC-NUE, when the number of nodes in the same area increases, it is more 
likely that among neighbouring nodes there are node where the processing cost is lower; 

• the higher the number of nodes in the same area, the higher the number of clusters formed, and 
therefore the bigger the amount of data that can be processed before they arrive to the Coordinator, 
reducing the energy cost. 

It may be inferred from the results that using the framework would be particularly energy conserving 
when data from different nodes have to be processed together, the processing is pretty complex, and the 
energy consumption or the initial energy is not uniform for the network. 

Figure 7. Percentage of energy conservation with respect to the ratio between processing 
cost and the cost to transmit 137 bytes of data, for cases UC-UE, NUC-UE and UC-NUE, 
for a density of 0.3 nodes per m 2 . Solid lines show energy conservation with respect to data 
processed only by the Coordinator; dashed lines show energy conservation with respect to 
data processed by every Cluster Head. 




0.5 0.7 0.9 1.1 1.3 1.5 
Processing cost to transmission cost ratio 



— ♦— C, case 1 
■ C, case 2 
* C, case 3 
-♦- CH, case 1 
-■- CH,case 2 
-A- CH,case 3 
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Table 4. Percentage values of energy conservation using the framework, for OpA and OpB, 
in cases UC-UE, NUC-UE and UC-NUE, comparisons C, CH and M, for a node density 
of 0.2 and 0.4 nodes/m 2 . 





Node density 
[nodes/m 2 ] 


Case 1 [%] 


Case2[%] 


Case 3 [%] 




C 


CH 


M 


C 


CH 


M 


C 


CH 


M 


OpA 


0.2 


19.5 


5.5 


11.7 


50.0 


29.3 


42.2 


47.5 


16.4 


33.6 


0.4 


25.6 


5.9 


16.5 


58.0 


35.3 


46.2 


56.7 


23.5 


43.7 


OpB 


0.2 


28.9 


0.1 


17.0 


33.3 


0.7 


19.9 


38.6 


1.6 


19.7 


0.4 


30.8 


0.1 


18.2 


37.8 


1.0 


21.4 


40.9 


3.5 


21.1 



In the considered framework we have not addressed the routing problem and we have assumed that 
the routing of the packets is working correctly. However, computation of the paths and re-computation 
in case of failure is for sure another cause of energy consumption. When the network experiences the 
failure of a node, that node must be bypassed and data addressed to it must be sent to the following node. 
To do it, an appropriate new routing path must be found. The number of packets exchanged to find a 
new routing path depends on the routing algorithm used by the network. If we consider a bad scenario 
where an average of 50 packets have to be sent among 10-hop distant nodes where each node is 2 m 
far from the other ones, in case of one node failure every hour, energy conservation would decrease by 
about 7.8%. It has been estimated that the decrease in energy conservation for OpA and for a density 
of 0.3 nodes/m 2 would be around 8.2% in case UC-UE, 12% in csae NUC-UE and 14.7% in case 
UC-NUE. 

Figure 8. Percentage of energy conservation for every area of the network for OpA, 
comparison C, in cases UC-UE, NUC-UE and UC-NUE, for a density of 0.3 nodes/m 2 , 
for an increasing distance of the area from the Coordinator. 
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A final observation may be made in regard to Figure 8, which depicts the percentage of energy 
conservation for OpA, comparison C, in cases UC-UE, NUC-UE and UC-NUE, for a density 
of 0.3 nodes/m 2 , in relation to the distance of each area from the Coordinator. As could be expected, the 
greater the distance of the sources from the Coordinator, the more energy conservation is derived from 
the use of the framework. 

Although they have not been reported, similar results have been obtained for all other cases. 

6. Conclusions 

In this paper we have studied the deployment of distributed applications in Wireless Sensor Networks 
and proposed a solution aimed at minimizing the impact of energy consumption on the network lifetime. 
We have considered a scenario where the nodes have dissimilar capacities in terms of sensing and 
processing. The resulting algorithm has been implemented to perform simulations in different scenarios 
and the results have been compared with alternative solutions. We observed significant improvements 
in terms of energy savings. We may therefore infer that using the framework described would be 
particularly energy conserving when the application encompasses the processing of data coming from 
different nodes, the processing is pretty complex, and the energy consumption of nodes as well as battery 
energy is not uniform over the entire network. 

References 

1. Yick, J.; Mukherjee, B.; Ghosal, D. Wireless sensor network survey. Comput. Networks 2008, 
52, 2292-2330. 

2. Atzori, L.; Iera, A.; Morabito, G. The internet of things: A survey. Comput. Networks 2010, 
54, 2787-2805. 

3. Singh, S.; Woo, M.; Raghavendra, C. Power-aware routing in mobile ad hoc networks. In 
Proceeding of the 4th ACM/IEEE International Conference on Mobile Computing and Networking, 
Dallas, TX, USA, July 1998. 

4. Chang, J.; Tassiulas, L. Energy conserving routing in wireless sensor networks. In Proceedigns 
of the Nineteenth Annual Joint Conference of the IEEE Computer and Communications Societies, 
INFOCOM 2000, Tel Aviv-Yafo, Israel, 26-30 March 2000; Volume 1, pp. 22-31. 

5. Park, J.; Sahni, S. An online heuristic for maximum lifetime routing in wireless sensor networks. 
IEEE Trans. Comput. 2006,55, 1048-1056. 

6. Liu, E; Tsui, C; Zhang, Y.J. Joint routing and sleep scheduling for lifetime maximization of 
wireless sensor networks. IEEE Trans. Wireless Commun. 2010, 9, 2258-2267. 

7. Ramanathan, R.; Hain, R. Topology control of multihop wireless networks using transmit power 
adjustment. In Proceedigns of the Nineteenth Annual Joint Conference of the IEEE Computer and 
Communications Societies, INFOCOM 2000, Tel Aviv-Yafo, Israel, 26-30 March 2000. 

8. Calhoun, B.; Daly, D.; Verma, N.; Finchelstein, D.; Wentzloff, D.; Wang, A.; Cho, S.; 
Chandrakasan, A. Design considerations for ultra- low energy wireless microsensor nodes. IEEE 
Trans. Comput. 2005, 54, 727-740. 



Sensors 2011, 11 



7419 



9. Ergen, S.C.; Varaiya, P. Optimal placement of relay nodes for energy efficiency in sensor networks. 
In Proceedings of the IEEE International Conference on Communications, ICC '06, Istanbul, 
Turkey, 1-15 June 2006. 

10. Tang, J.; Hao, B.; Sen, A. Relay node placement in large scale wireless sensor networks. Comput. 
Commun. 2006, 29, 490-501. 

11. Gao, Q.; Blow, K.; Holding, D.; Iqbal, Q.R. Radio range adjustment for energy efficient wireless 
sensor networks. In Proceedings of ' SpringSim' 07 , Norfolk, VA, USA, 25-29 March 2007. 

12. Howitt, I.; Wang, J. Energy balanced chain in wireless sensor networks. In Proceedings of the 
IEEE Wireless Communications and Networking Conference, WCNC 2004, Atlanta, GA, USA, 
21-25 March 2004. 

13. Heinzelman, W.; Chandrakasan, A.; Balakrishnan, H. An application-specific protocol architecture 
for wireless microsensor networks. IEEE Trans. Wireless Commun. 2002, 1, 660-670. 

14. Hong, T; Wu, C. An improved weighted clustering algorithm for determination of application 
nodes in heterogeneous sensor networks. J. Inf. Hiding Multimedia Signal Process. 2011, 
2, 173-184. 

15. Spiess, R; Karnouskos, S.; Guinard, D.; Savio, D.; Baecker, O.; de Souza, L.M.; Trifa, V. 
SOA-Based integration of the internet of things in enterprise services. In Proceedings of the 
International Workshop on Cloud Services, IWCS 2009, Los Angeles, CA, USA, 10 July 2009. 

16. Wang, Q.; Hempstead, M.; Yang, W. A realistic power consumption model for wireless sensor 
network devices. In Proceedings of the Third Annual IEEE Communications Society Conference 
on Sensor, Mesh and Ad Hoc Communications and Networks, Reston, VA, USA, 25-28 Sepember 
2006. 

17. ILOG. CPLEX User's Manual. Available online: http:// yalma.fime.uanl.mx/cplexl 1 -manual/ 
Content/Optimization/Documentation/CPLEX/_pubskel/ XPlatform/User _man250.html (accessed 
on 9 June 2011). 

18. FICO. Xpress Optimization Suite. Available online: http://www.fico.com/en/Products/DMTools/ 
Pages/FICO-Xpress-Optimization-Suite.aspx (accessed on 9 June 2011). 

19. Chipcon. SmartRF CC2420 Datasheet. 2.4 GHz IEEE 802.15.4/ZigBee-Ready RF Transceiver. 
Available online: http://inst.eecs.berkeley.edu/ csl50/Documents/CC2420.pdf (accessed on 9 June 
2011). 

20. IEEE. Std 802.15.4™, 2006. Wireless Medium Access Control (MAC) and Physical Layer 
(PHY) Specifications for Low-Rate Wireless Personal Area Networks (WPANs). Available online: 
http://www.google.com.hk/url?sa=t&source=web&cd=l&ved=0CB0QFjAA&url=http%3A%2F% 
2Fprofsite.um.ac.ir%2F hyaghmae%2FACN%2FWSNMACl.pdf (accessed on 9 June 201 1). 

© 2011 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article 
distributed under the terms and conditions of the Creative Commons Attribution license 
(http://creativecommons.Org/licenses/by/3.0/). 



